SRB STATES AND NONEQUILIBRIUM STATISTICAL 
MECHANICS CLOSE TO EQUILIBRIUM. 



by Giovanni Gallavotti* and David Ruelle** 

Abstract. Nonequilibrium statistical mechanics close to equilibrium is studied using 
SRB states and a formula [10] for their derivatives with respect to parameters. We write 
general expressions for the thermodynamic fluxes (or currents) and the transport coeffi- 
cients, generalizing the results of [4], [5]. In this framework we give a general proof of the 
Onsager reciprocity relations. 

There is currently a strong revival of nonequilibrium statistical mechanics, based on 
taking seriously the nonlinear (chaotic) microscopic dynamics. One natural idea in this 
direction is to use nonequilibrium steady states, which are defined to be the SRB states 
for the dynamics (see below). 

The idea of using SRB states eventually led to useful results only recently, when it 
was combined in [6] with reversibility of the dynamics to obtain a nontrivial fluctuation 
formula for the entropy production. In [5] the same ideas were applied to prove the Onsager 
reciprocity relations, and the fluctuation-dissipation formula for a rather special class of 
models. The analysis dealt with examples rather than the general situation, and relied on 
an unproven conjecture on Anosov systems. In this note we generalize [5] and, using [10], 
give results that can be proved rigorously for Axiom A diffeomorphisms (this is a strong 
chaoticity assumption, see [12], here we skip technical details). 

To be discussed is the classical microscopic description of the time evolution of a 
physical system. (No large system assumption will be made) . We let the accessible phase 
space be a compact manifold M. The time evolution is given by iterates of a diffeomorphism 
/ of M (discrete time case) or by integrating a vector field F on M (continuous time case). 
We choose a Riemann metric on M, and let dx be the corresponding volume element. 

In nonequilibrium statistical mechanics, the time evolution typically does not preserve 
any measure which has a density with respect to dx. Let m(dx) = m(x)dx be a probability 
measure (with density to), and f* k m the direct image of m by f* k ; any weak limit for 
n — > oo of 
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is an /-invariant probability measure pf on M. If furthermore pf is ergodic we may say 
that it is a natural nonequilibrium state (the SRB states are special examples of this, [3]). 

An infinitesimal change 8f of / corresponds to an (infinitesimal) vector field X = 
Sf o f~ x on M and an easy formal calculation gives 

oo 

%(<&) = $>/<grad(<I>o /"),*> (1) 

n=0 

where $ is a smooth test function and (■ , • ) is the scalar product of a tangent vector and a 
cotangent vector to M. Equation (1) expresses the change in the expectation value of the 
observable $ when the system is subjected to a force X. Replacing X by a time dependent 
force Xt we have for the change of expectation value of the observable $ at time s the 
analogous formula 

8 sPf (<S>) = £>/<grad(* o f-%X t ) (2) 

t<s 

Note that the condition t < s corresponds to the "causality principle" and (2) can 
be used to derive Kramers-Kronig dispersion relations. For Axiom A diffeomorphisms, (1) 
and (2) can be proved rigorously*: the left-hand side is a derivative, and the right-hand 
side a convergent series. 

In the case of continuous time systems described by a differential equation x = F(x) 
and by the corresponding flow (/*), an infinitesimal variation SF = X generates a variation 
in the expectation value Pf(®) (for the natural nonequilibrium state) given by: 

r+oo p 

5p F ($) = dt J p F (dx)(gmd x (<f>of),X(x)) (V) 

(A rigorous proof for Axiom A flows has not been given yet). 

The entropy production associated with the diffeomorphism / is defined by** 

e f = Pf( (T f) . (T / = - lo g J / ( 3 ) 

where J/ is the absolute value of the Jacobian of / with respect to the Riemann volume 
element dx. In the continuous time systems, we let 

eF = Pf{o~f) , o~f = — divF (3') 

(see [5]). 



* See [10]; the proof given in [10] for the time dependent case assumes that the pertu- 
bation has finite support in time. 



For a discussion of entropy production see [9] . 
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From now on we shall fix / such that Pf(dx) = dx, i.e., pf is the Riemann volume 
element (if pf has smooth density this can be achieved by a change of metric). Note that 
in particular this implies that (1) can be rewritten as 

00 

*P/(*) = j>/(*o^(-divX)) (4) 

n=0 

In this nondissipative situation we have log Jf = 0, hence e/ = 0. If we write as before 
X = 5f o / _1 we obtain to second order in X, using (3) and (1), 

e f+Sf = (Pf +Spf)(cTf+Sf) 

= ^p f ((dwX) 2 )-Sp f (divXof) 

1 00 
= - W ((divl) 2 ) - ^ P/ (grad((divX) o n,X) 

1 00 

= -p f ((dwX) 2 ) + PfddivX) o r. divX) (5) 

Z n=l 

- Pf((divX)or.dwX) 
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The same analysis leads, in the continuous time case (with 5F = X, and pp(dx) = dx) 

to: 

^ r+00 

e F+SF = g J dt Pf(°f /* ' <*f) 

1 f 

- dt dxdivXifx) ■ divX(x) (5') 

2 j-00 7 
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see [4], [5]. 



We shall now relate the above expressions for the entropy production to the definition 
of the thermodynamic forces X a , and the conjugated thermodynamic fluxes J a as they 
appear for instance in [8]. We begin by an informal discussion, and assume, as usual in 
applications, that / or F depends on parameters E a , so that we may write (to first order) 

x = Y J y^E a 

OL 

We identify the thermodynamic forces X a with the parameters E a . Considering first the 
continuous time case, we follow [4], [5], and we define the thermodynamic flux conjugated 
to E a as 

, d . 
Jo. = Pf+sf{wttO-f+5f) 



Since da/dE a is a divergence, p(da/dE a ) = and we have 

Jo. = 5p F (-divV a ) + h.o. 



= J dtj PF (dx)(gvad x ((-dwV a ) o f),X(x)} + h.o. 

From now on we neglect higher order terms and (using integration by parts, since pp(dx) = 
dx) we write 

dt J p F (dx) (div x X)(div / i x V Q ) 

In the discrete time case we define the thermodynamic flux only to leading order in 
the E a by 

oo 

J a = - Pf ((dwV a )(dwX)) + Pf((divV a ) o /». (divX)) 

n=l 

With these definitions J a depends only on the application of X in the past (causality, 
cf (2)) and the entropy production (to second order) is 

ef+sf or eF+SF = ^ XqJa 

a 

These conditions uniquely determine the J a . Notice that the formulae for J a involve only 
the divergences of X and V a . 

To continue the discussion, we assume that there is a (sufficiently large) Banach space 
B of functions $ : M — > R such that 

p/($)=0 if $G6 

and for some constant C 

J2\Pf(^ ° f k -®)\ < C||^IIb||*||b if $,*e£ (6) 

kez 

[This is the discrete time case, the continuous time case is similar. If / is an Anosov 
diffeomorphism* we can take for B a space of Holder continuous functions on M. Similarly 
for Anosov flows]. 

From now on, we assume that divX is in the Banach space B just introduced, and we 
may define X e B and J E B* (the dual of B) as follows: 

X = -divX G B (7) 

1 oo 

(J, *) = 7jP/("divX $) + ^ P/ ((-divX) o /». $) (8) 

Z n=l 



* An Axiom A diffeomorphism / preserving dx is an Anosov diffeomorphism. 
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for discrete time, and 

(J,Q)= ^ jf p f ((-divX)o/ ( .$) (8') 

for continuous time, where $ G £>, and (•,•) is the pairing £>* x B — > C. Note that if 
divV a G £> we have 

Ja = (J, -divK) 

and that the entropy production is 

ef+sf or e F +8F = (J, X) 



With the above notation and assumptions we may write J = LX , where L is, in view 
of (6), a continuous linear map £>-+£>*. If we define a unitary operator U on L 2 (pf) by 
= $ o / and write f7 = e* a dP(a) we have 

J] e- ika p f {& o f k . $) = 2tt-^-($, rfP(a)$) L2 

so that the quadratic form associated with L satisfies 

(Ld>, $) = i £ P/ (<& o $) = tt dP( a )$) L2 | a=o > 

fcez " (! 

In particular, this quadratic form is > 0. 

The formulae obtained up to now hold under the only assumption of closeness to 
equilibrium. If we make the further assumption of (microscopic) reversibility, we shall 
obtain a symmetry property of L called Onsager reciprocity. For simplicity we discuss 
only the discrete time case. 

We say that the dynamics is reversible if there exists a diffeomorphism i : M — > M 
such that i 2 = identity, iof = f~ 1 oi. We have then also i*pf = pf. [Note that pf, i.e., the 
Riemann volume, is mixing by (6), hence /-ergodic. Since i*Pf is absolutely continuous 
with respect to pf and satisfies f*(i*pf) = i* (/ _1 )* Pf = (i*Pf) we have i*pf = pf by 
ergodicity] . Assuming reversibility we may define e$ = $ o j for $ G £>, and we find 

OO OO 
n=l n=l 

OO 

= - P/ (e*. e¥) + £ p/(e$ o /"». e¥) = (L(e<&), e¥) 

Z n=l 

The relation 

(Ltf,$) = (L(e$),e#) 
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is a form of the Onsager reciprocity relation as we shall see in a moment. Note that 
reversibility was assumed only for / (i.e., at equilibrium), the perturbation 5f is arbitrary. 

To obtain a more familiar form of the entropy production formula (see [8]), we assume 
that B has a basis (<E> a ) with a corresponding system (0 a ) of elements of B* such that 
(0 a , $0) = <W (see [11]). Write 

X a = (0 Q , X) = (0 Q , -divX) 
00 

Ja = (J, $a) = -p/("divX, $ Q ) + ]T P/(("divX) O /». d> a ) 

In particular 



n=0 



and 



-dM = ]T(0 Q , -divX). $ Q = J2 X «- ^ 

a a 

e f+5f = (J,X) = E ^ Q 
a 

00 

j° = 2 E ^m^) + E E /"• 

/3 n=l /3 

To avoid convergence problems suppose that finitely many Xp only are nonzero. Then 

Ja = ^ L a pXp 

P 

where 

1 

L *P = -Pf($f3$a) + Pf@0 ° 

n=l 

Suppose again reversibility of the dynamics, and let the basis (<E» a ) of B be such that 
e$ Q = $a o / = e a $a with e a = ±1. Then 

Lap = tatpLfia 

which is the usual form of the Onsager reciprocity relation. 

We conclude by sketching an example, see [1], of the formalism described above. Let £ 
be a surface of constant negative curvature, and genus g, with the g — 1 automorphic forms 
(j) a (z)dz. We normalize these forms so that they are ortho normal in the space L2(T£) in 
the natural scalar product, [7]. 

We can then consider the hamiltonian equations of the motion of a particle on £ 
subject to the external force generated by the "electric" field 8 such that £ x dx + £ y dy = 
Re Y. E oi¥a{z)dz. 
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We also impose, via Gauss' principle, [5] , that there is a thermostat force that keeps 
the kinetic energy constant (and equal to 1/2) in spite of the field's action. Thus the 
"only" efffect that the fields will have on the flow is that currents flowing "around" the g 
"holes" of the surface will be established. But the flow being a mixing Anosov flow on a 
compact surface (because of the gaussian constraint) it will result that a stationary state 
will be reached and the latter will be the SRB distribution, [2] . 

The equations of motion can be written explicitly and one finds in particular that the 
entropy creation rate at the point q, p of the phase space is 

The transport coefficient can also be explicitly computed, and one finds L a p = |5 a /3, see 
[!]• 



References. 

[I] F.Bonetto, G. Gentile, V.Mastropietro. "Electric fields on a surface of constant 
negative curvature." Preprint, 1996. 

[2] R.Bowen, D.Ruelle. "Ergodic theory of Axiom A flows." Inventiones Math. 
29,181-202(1975). 

[3] J.-P.Eckmann and D.Ruelle. "Ergodic theory of chaos and strange attractors" Rev. 
Mod. Phys. 57,617-656(1985). 

[4] G.Gallavotti: "Extension of Onsager's reciprocity to large fields and the chaotic 
hypothesis." Phys. Rev. Letters. 77,4334-4337(1996). 

[5] G.Gallavotti. " Chaotic hypothesis: Onsager reciprocity and fluctuation-dissipation 
theorem." J. Statist. Phys. 84, 899-926(1996) 

[6] G.Gallavotti and E.G. D.Cohen. "Dynamical ensembles in nonequilibrium sta- 
tistical mechanics." Phys. Rev. Letters 74,2694-2697(1995). "Dynamical ensembles in 
stationary states." J. Statist. Phys. 80,931-970(1995). 

[7] I.M.Gel'fand, M.I.Graev. I.I.Pyateckii-Shapiro. Representation theory and auto- 
morphic functions. Saunders, Philadelphia, 1969. 

[8] S.R.de Groot and P.Mazur. N on- equilibrium statistical thermodynamics. Dover, 
New York, 1984. 

[9] D.Ruelle. "Positivity of entropy production in nonequilibrium statistical mechan- 
ics." J. Statist. Physics, 85,1-25(1996). 

[10] D.Ruelle. "Differentiation of SRB states." To be published; the preprint is 
archived inmp_arc@math.utexas.edu, 96, #499. 

[II] I. Singer. Bases in Banach spaces I. Springer, Berlin, 1970. 

[12] S.Smale. " Differentiable dynamical systems." Bull. Amer. Math. Soc. 73,747- 
817(1967). 



7 



